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Abstract 

We study the theory of non-relativistic matter coupled to the non-Abelian U{2) 
Chern-Simons gauge field in (2+1) dimensions. We adopt the mean field approximation 
in the current- algebra formulation already applied to the Abelian anyons. We first show 
that this method is able to describe both "boson-based" and "fermion-based" anyons 
and yields consistent results over the whole range of fractional statistics. In the non- 
Abelian theory, we find a superfluid (and superconductive) phase, which is smoothly 
connected with the Abelian superfluid phase originally discovered by Laughlin. The 
characteristic massless excitation is the Goldstone particle of the specific mechanism of 
spontaneous symmetry breaking. An additional massive mode is found by diagonalizing 
the non-local, non-Abelian Hamiltonian in the radial gauge. 
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1 Introduction 



The dynamics of anyons - particle with fractional statistics in (2 + 1) dimensions |jj - has 
been considerably investigated in the past few years. These collective excitations can arise 
in planar condensed-matter systems like the fractional quantum Hall effect and the high- 
temperature superconductivity [Qj. An effective field theory for anyons is obtained by cou- 
pling non-relativistic matter particles - either bosons or fermions - to the Abelian Chern- 
Simons gauge field, which provides the statistical interaction [|J ||. 

The remarkable property of superfluidity is exhibited by anyons in the thermodynamic 
limit at constant density |J : the anyon fluid possesses a ground state with uniform density 
and a massless longitudinal excitation. This is a Goldstone mode which gives rise to super- 
conductivity by the usual Higgs mechanism when coupled to the physical electromagnetic 
field. This theory was originally proposed by Laughlin |§ for explaining the high-temperature 
superconductivity of cuprates || . However, the explicit breaking of P and T symmetries by 
the fractional statistics [|IJ has not been confirmed by the experiments so far |7]]. 

Independently of its physical application to high-temperature superconductivity, we be- 
lieve that the anyon superfluid is very interesting and deserves a deeper analysis. Few 
non-perturbative, semiclassical, ground states are known in field theory, thus any new one 
is worth understanding for its own sake. This may find wider applications than the origi- 
nal physical problem, as it has occurred to spontaneous symmetry breaking. Actually, the 
Anyon fluid is closely related to the usual superfluid, because both exhibit the spontaneous 
breaking of the U(l) global symmetry related to particle number conservation. 

In this paper, we show that the anyon superfluidity also arises in the presence of a 
non- Abelian Chern-Simons interactions. We consider the simplest case of particles having a 
isospin 1/2 quantum number with U(2) gauge symmetry. The application of this theory to 
some (2 + l)-dimensional physical problems has been discussed in ref. ||. Here, we solve it 
in the mean field approximation, describe the specific mechanism of spontaneous symmetry 
breaking and study the low-energy excitations. 

In section two, we review the mean field approximation || which allows to describe the 
non-perturbative ground state of the Abelian anyons. One assumes self-consistently that 
the matter density is spatially uniform and obtains a uniform magnetic field by the Chern- 
Simons Gauss law, (p) = —(B)/n , where k is the coupling constant. Thus, the particles 
uniformly fill up the Landau levels determined by this mean magnetic field. The quadratic 
fluctuations around the mean field can be described by using the Dashen-Sharp current 
algebra formalism [10]; their diagonalization by a Bogoliubov transformation produces a 
relativistic longitudinal excitation at low energy, as in the familiar case of the superfluid 
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This is by far the simplest method for describing anyon superfluidity |J. 

While reviewing this method |J, we clarify one property of the Chern-Simons interaction 
in the Hamiltonian formulation. This long-range, topological interaction produces non- 
trivial boundary effects, whose strength depends on the type of mean field ground state. 
These boundary effects can be removed by normal-ordering the Hamiltonian, but produce 
an effective local interaction, which is ground-state dependent. This property is crucial 
for describing both "boson-based" and "fermion-based" anyons within the current algebra 
approach. Actually, the two descriptions of anyons are more accurate for fractional statistics 
9/tt ~ and 9/ir ~ 1, respectively, and agree at the mid-point of semions (0/tt = 1/2): their 
combination yields a consistent approximation for all values of the statistics. In particular, 
we obtain the approximate second-order ground-state energy. 

In section three, we extend this method to the U(2) non-Abelian Chern-Simons interac- 
tion with two independent coupling constants, Ku(i) = « and ksu(2) = k. The mean field 
approximation produces two copies of Landau levels for isospin-up and isospin-down matter, 
which have opposite contributions to the mean iso-magnetic field. For kR > 0, the ground- 
state configuration corresponds to equal populations of spin-up and spin-down particles, and 
to a vanishing iso-magnetic field; if 1/k — > 0, this ground state is P and T invariant because 
parity-violating effects cancel between the two populations - only excitations can break P 
and T explicitly. Models of this kind have been discussed in refs. flT2||. 



Here we describe a different phase of the system, which exists for kk < and < 
4/ 1 «; |, and has a ground state with maximally unbalanced populations. This phase is con- 
tinuously connected to the Abelian theory by tuning 1/|k| — > 0. The ground state breaks 
spontaneously the U{2) global symmetry to a U{1) subgroup, as in the Standard Model of 
electroweak interactions We find it interesting that the low-energy dynamics of a non- 
Abelian gauge theory can be solved in closed form in a toy model for spontaneous symmetry 
breaking. Moreover, the dynamics of non-Abelian anyons has not been much investigated so 
far ill. 

In section four, we discuss the low-energy collective excitations above the mean-field 
ground state. The quadratic expansion of the Hamiltonian, written in terms of non-Abelian 
currents, consists of two independent parts, corresponding to matter-density and isospin- 
density fluctuations, respectively. The former fluctuations are massless and similar to those 
of the Abelian anyon fluid. The latter have a non-trivial, yet solvable, non-Abelian dynamics; 
we solve explicitly the Gauss law constraint by using the radial gauge Y%=\ x l A^ = 0, which 



maintains manifest rotation invariance and breaks translation invariance 15 . The non- 



Abelian Hamiltonian resembles a non-local deformation of the Landau-level Hamiltonian, 
because isospin-flip excitations feel the mean iso-magnetic field and, moreover, self-interact. 
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We obtain the complete spectrum and show that this is gapful and discrete, with Gaussian 
fall-off of correlations. Special care is paid to the gauge invariance of the set of physical 
states, which is actually translational invariant. This massive excitation does not spoil the 
Laughlin superconductivity mechanism, because both low energy excitations become gapful 
upon coupling to the physical electro-magnetic field. 

In section five, the Abelian and non-Abelian anyon superfluidities are explained in terms 
of the spontaneous breaking of the global gauge symmetries, and the specific mechanisms are 
compared with those of the Higgs and Standard Models of four dimensional gauge theories 
|13| . Finally, in the conclusion, we discuss other possible physical applications of the non- 
Abelian anyon fluid. In the appendix, we collect some additional informations on the eigen- 
functions of the non-Abelian Hamiltonian. 



2 The mean field approximation in terms of currents 
2.1 Hamiltonian and Abelian current algebra 

In this section, we review the mean field for the Abelian anyon fluid in the current algebra 
approach of ref . [[J . The Lagrangian for non-relativistic matter coupled to the Abelian Chern- 
Simons gauge field is ||, 

C = itfDo* - ^(A*) f A* + ^A a d^ , (2.1) 

where A = <9 M + iA^ is the covariant derivative of the gauge field and \I/ is the non- 
relativistic matter field[|. The equation of motion for the matter field can be used to derive 
the Hamiltonian, 

n = J rf 2 x -L(A*) f (A*) , (2.2) 

while the equations for the gauge field are, 

— F12 = B = eijdiAj = p (Gauss law), 

Kj 

F 0l = d Ai + diA = -- ey J j . (2.3) 

K 

The conserved, gauge- invariant matter current J M = (p, J 1 = Jj) is given by 

p = tfttf j i = J_(^t D .^ _ (£>.tf)ttf) . (2.4) 

Lim 

* We choose units c — h — 1 and set the electric charge equal to one. We use the metric rj^ = 
diag(l, — 1, — 1), while for two-dimensional expressions we use 5{j, i.e. — (Aq,— A 1 — —AA and x 11 — 

(x , X — Xi ) . 
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The Chern-Simons field has no local physical degrees of freedom, thus it can be solved in 
terms of the matter field at equal time by choosing a complete gauge fixing. The Coulomb 
gauge diAi = has been often used for non-relativistic theories; here, we prefer the (spatial) 
radial gauge [ |i~5|l , 

x% = , (2.5) 

because it is better suited for the non-Abelian theory discussed in the next section. Actually, 
any choice of gauge is equivalent for the Abelian theory, because it will be described in terms 
of gauge-invariant quantities. The solutions of gauge field equations (|2.3|) in the radial gauge 



are 



A 
A, 



,IJ_ e .,7. 

1 j 

eijX J p 



(2.6) 



k 1 + x ■ d ~ %r 

and do not actually involve time derivatives. A precise meaning of the operator (x ■ d) 1 
is not important here, and will be discussed in section four. By using (|2.6| ), we can write 
the Hamiltonian in terms of matter fields only, and quantize it by requiring the bosonic 
commutation relations 

*(x,t),*t(y jt )] =5(x-y) . (2.7) 

Let us now choose the current Jj and the density p as basic variables. Their algebra can 
be computed by using (|2.7| ) and (|2lf), and reads, 



b(x),p(y)] 
[p(x),Ji(y)] 

[^(x),J,(y)] 



o, 

1 d 



im dx 
1 



(<5(x-y)p(x)) 



im 



A (5(x _y )j4(x)) _|_ (5(x _y )Jj(y)) 



{21 



An important property of this algebra is its independence of the Chern-Simons coupling 
constant, which only appears in the representation of the algebra (the states) and in the 
normal- ordered Hamiltonian [[J. The Hamiltonian (|2.2|) can also be written in terms of 
currents as follows [flOfl , 



r 1 1 

H = / d 2 x (dip + 2imJi)^ - (dip + 2imJi) . 

J 8m p 



(2.9) 



2.2 Mean field approximation 

Let us assume that the ground state \Q) has a spatially uniform density, 

(fi|p(x)|ft) =po, (ft|Ji(x)|ft) = , 



(2.10) 
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which corresponds to a uniform magnetic field, 

(n\B(x)\{i) = B = --p , (2.11; 



K 



by the Gauss law (|2.3| ). Next, we seek for self-consistency of this hypothesis in the approx- 
imate quantum theory. By decoupling matter (p) and field (B) fluctuations in Q2.ll 
quantize the field $ in the external average magnetic field Bq\ 



we 



H - w(o) = j d 2 x ^-(Dl°h)\D^) , (2.12) 

where Z)j = d{ — iAf\ and A[ = e^x-* po/2n. This is the well-known Hamiltonian of the 

Landau levels in the so-called symmetric gauge [fL6|| . The one-particle energy of the n-th 
Landau level is, 

e n = ^(n + h, (2.13) 
m 2 

and the eigen-functions of the lowest level are 

<Mx) = * (£)' e-l^ 2 /^ 2 , (W + ** 2 ) , (2.14) 



where A = y2/eB is the magnetic length, and £ is the angular momentum. Note that these 
eigen-functions satisfy 

(£)} 0) +i£)f)^ = 0. (2.15) 

The angular momentum orbitals ( |2.14| ) have degenerate energy, and their number is 
(BoA/2n) in a finite domain of area A (independent of n). 

The mean field hypothesis is self-consistent for all the ground states of Ti.(°> with N 
particles which have uniform density. For bosonic matter, these have been found in ref.0, 
and correspond to filling each Landau orbital of the lowest level with the same number n of 
particles. Actually, using ( |2.14jj2.15| ), one can compute 



lim (fi|p(x)|fi) = — = po , lim (Q\Jihc)\Q) = . (2.16) 



These values for p agree with the Gauss law ( |2.3| ), provided that k = —n/2Tr. Therefore, the 
mean field approximation is self-consistent for these integer values of the coupling constant. 
The ground-state energies, for a system of area A, are obtained from ( |2.13| ), 

= to N = po A , (boson — based anyons, k = ) . (2-17) 

nm V 2n J 

The next order of the mean-field approximation is given by the quadratic fluctuations 
of the density and the current. The Hamiltonian ( O ) must be expanded quadratically and 
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normal-ordered in the thermodynamic limit N — > oo. The latter limit involves some subtle 
boundary effects which actually determine the strength of the effective local interaction of 
fluctuations. Actually, a more precise expression of ( |2.16| ) for large, but finite, N can be 
obtained from (gg , glg ) |7|, 



(ft|p(x)|fi)~p e 



'N 



— 7T X 



(fi|Ji(x)|fi) 



2m 



e l0 a,-<fl|/>(x)|fl) . (2.18) 



Namely, the density of a filled, finite, Landau level has the shape of a droplet, with a chiral 



edge current, due to eq. (|2.15|) . The contribution of this edge current to the ground-state 
value of the Hamiltonian in the form ( p.9[ ) is non- vanishing for N — > oo, and correctly gives 
the ground-state energy ( |2.17| ). This boundary effect can be removed by rewriting the 
Hamiltonian. Using an algebraic identity of the Bogomol'nyi type H, 



ie ij (Di^y Dfl + me ij diJ j + Bp = 



the Hamiltonian (|2.2|) can be rewritten 



n 



d 2 x 



1 

2m 



(A*) f (A*) + i« ^-(A^A* - 



a 



(2.19) 



(2.20) 



for any value of a. For a — 1, the derivative terms in (12.20] ) vanish on the ground state, 
due to eq. ( |2.15| ); thus, there are no boundary effects. The ground-state energy ( |2.17| ) is 
given by the local term p 2 only. The new expression ( 2.20Q of the Hamiltonian can be easily 
normal ordered in the thermodynamic limit as follows: 



:H:=n- (Q\H {0) \Q) =H-— I d 2 *pl 

nm J 



(2.21) 



Therefore, the Hamiltonian (|2.20j) with a — 1 is adopted for studying the quadratic fluctua- 
tions. Note that the normal-ordering procedure has produced an effective local interaction, 
which is the mean field approximation of the long-range "statistical repulsion" of anyons. 
This statistical repulsion generates a positive energy density Eq > as in the case of free 
fermions. 

Finally, we note that an attractive local interaction (—gp 2 /2) can also be included in the 
Hamiltonian for the anyon fluid ( |2.2| ), ( |2.20| ) Pfl ||18j| . The previous analysis can be extended 
for generic values of g < l/(m|/t|) . At the "self-dual" point g = l/(m|/c|), the local attraction 
exactly balances the statistical repulsion, and there is phase transition for the anyon fluid: 
non-trivial classical solutions with i^o = were found in ref . H and conformal invariance was 



shown to hold to three- loop order in ref. [18]. The nature of the other phase g > l/(m|/c|) is 
not presently understood. 
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2.3 Quadratic fluctuations 



Following the approach of ref. ||, we study the quadratic fluctuations using the variables 
(p, Jj), satisfying the algebra (|2.8|) . Actually, we are only interested in representing this 
algebra to leading order in the fluctuations, as well as expanding the Hamiltonian ( 2.20Q to 
quadratic order. To this effect, we introduce a small parameter e which keeps track of the 
size of fluctuations, 

p(x) =p + ep(x) + 0(e 2 ) , J(x) = eJ(x) + 0(e 2 ) . (2.22) 

By inserting this expansion in the second of the current commutators ( [2.8| ), we obtain 

Po f d 



P(x),j(y) 



(2.23) 



This shows that e 2 is of order 0(h), so that we can neglect the fluctuations in the r.h.s. of 
the commutators. The third commutator in (|2.8|) can be similarly estimated: 



J<(x), Jj{y) 



0(e 



. 



(2.24) 



In general, the use of this e-expansion yields consistent results for multiple commutators of 
the approximate algebra. 

This approximate algebra, given by (|2.23|) , ( j2.24j ) and [p(x), p(y)] = 0, can be represented 
by a bosonic canonical field <ft, satisfying 0(x), </^(y) = 5(x — y), as follows, 



P(X) = y/to(<f> + 0^) 



Ji(x) 



2im 



0,( 



(2.25) 



(Note that the zero mode of <fi is absent due to the condition J c/ 2 x p(x) = N. ) The 
Hamiltonian ( p.20|) , with a = 1, can be written in terms of currents, similarly to ( |2.9|) , and 
then expanded to quadratic order in the fluctuations. The result is M 



rf 2 x: 



2m 



-j—{K\ + tK 2 )\Ki + iK 2 ) + — p 2 (x) 
4p n 



where Ki = d^p + 2imJi. By inserting the Fourier modes 

d 2 p 



2tt 



dp, dq 



2 (p-q) 



one obtains, 



W (2) 



c/ 2 p : Apapdp + B(a p a- 

27r P0 n TTPO 



p T" ftpGl_p, 



2m nm 



nm 



(2.26) 



(2.27) 



(2.28) 
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The terms coming from the repulsive interaction can be normal-ordered by the Bogoliubov 
transformation 



cosh x a P + sinh \ «- P , X = x(p) , 
aL p = sinhx a p + coshx aL p . (2.29) 

This is an 5*0(1, 1) rotation in the (a p , aL p ) space which preserves the commutation relations. 
The result is 

H {2) = AE^ + J d 2 p E p apa p , (2.30) 

with 



|p| /27rpo V 2 p^q | , ( qi x 

£ p = — W + — — ► v a p , 2.31 

m V n 4 



and 



= ~(27rP 7 4m ( ^"V T 2 / ' ~ = ~7T ^ 



(27r) 2 J Am \ p 2 y p : 
The quadratic fluctuations show a massless longitudinal excitation with sound velocity 



„. = I, *ZH , boson - based anyone « = - » i = I) (2.33) 

The use of the Bogoliubov transformation || makes manifest the striking similarity be- 
tween the anyon fluid and the usual superfluid JTTJ] . The latter is the canonical example 



for spontaneous breaking of a global symmetry , the U(l) symmetry for particle number 
conservation. Actually, the same spontaneous symmetry breaking occurs in the anyon fluid, 
because the Bogoliubov rotated ground state, satisfying a p |f2) = 0, does not have a well 
defined particle number. The broken symmetry is the global U(l) subgroup of the gauge 
groupQ. Therefore, the anyon fluid gives an interesting new realization of the Goldstone 
mechanism in non-relativistic field theory. Note that the "microscopic" mechanism leading 
to (p) = po is different from the Bose-Einstein condensation, and that there is no Higgs phe- 
nomenon associated to the Chern-Simons field. We shall discuss these differences in section 
five, together with the results of the non-Abelian case. The anyon fluid becomes a supercon- 
ductor when is coupled to an external electro-magnetic field, because the massless mode 
gives mass to the photon by the usual Higgs mechanism. 



2.4 Fermion- based anyons 

It is interesting to extend these results to fermion-based anyons. Suppose now that the 
matter field \I/ satisfies canonical anti-commutation relations. The mean field is again self- 
consistent (eq. ( p.llD ) for uniform fillings of the Landau levels, because each filled level 
' Besides, it is the unique part of the gauge symmetry which can break [19]. 
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contributes a constant value to the density, away from the boundary [P0| . Due to Fermi 
statistics, we can put two spin-1/2 fermions per Landau orbital, at most, and uniformly fill 
the lowest n/2 Landau levels, where n = 2p + a, a — 0, 1; if n is odd (er = 1), we fill the 
top level with one electron per orbital. The resulting ground-state density is again given 
by (|2.16| ), and the allowed values of the Chern-Simons coupling constant are K = —n/2-K , 
which correspond now to the fractional statistics 9/ir — 1 — 1/n. The ground-state energy, 
obtained by (|2.13|) , is: 

^ k=0 I -T&B- ( 1 + r?r I , n odd, 







2m ' 








2m V 


+ n 2 




n 


K = 









fermion — based anyons, k = ) . (2.34) 

This ground-state energy oscillates between even and odd values of n and correctly repro- 
duces the energy of the filled Fermi sea for k — » oo. 

The Hamiltonian must be normal-ordered differently from ( [2.21| ), because the ground- 
state energy is higher for fermion-based anyons than for boson-based ones. Again, we can 
dispose of the boundary terms in the ground-state expectation value of H by choosing the 
parameter a in eq. ( [2.20D which gives vanishing derivative terms. This is found to be 
a = n/2 + (a/2n) by using some equations similar to ( |2.15| ). As a consequence, fermion- 



based anyons have an effective local repulsion different from the boson-based ones. The 
discussion of quadratic fluctuations is the same as in the previous bosonic case, because the 
current algebra is independent of the statistics. We obtain the Hamiltonian ( |2.28| ) with 



' 2m m \ n z J 2m \ n A 

leading to a massless mode with sound velocity, 



^ y/ptfrr , n even 



n 



v _ = < , / ; TV fermion — based anyons , k = , — = 1 



(2.36) 

This value is different (always lower) than the result of ref.|J] for fermion-based anyons, 
because these authors considered spinless fermions. In the free fermion limit, v s — > Vpj\[2, 
where vp is the velocity of particle-hole excitations at the Fermi surface: thus, the mean 
field approximation picks up one particular value of the continuum of massless particle-hole 
excitations with velocities < v s < vf- Note also that this approach gives the same result for 
boson-based and fermion-based anyons at the common midpoint of semions, with statistics 
9/n = 1/2, for both and v s (eqs. (|2~T7|J235| ) and respectively). This 

shows that the current-algebra approach can describe both types of anyon constructions, 



9 



with greater accuracy in the regions 9/tt ~ and 9 /tt ~ 1, respectively, corresponding to 
1/ k, = 2-K/n — > in both cases 0. 



2.5 Ground- state energy 



The best approximation for the ground-state energy + AE^ is obtained by combining 
the boson-based expression for < 9/n < 1/2 and the fermion-based one for 1/2 < 9/n < 1. 
As discussed in ref. ||, the second-order contribution AE^ ( |2.32| ) is ultraviolet divergent 
and must be regularized by allowing a finite size to anyons, a = 1/A, where A is the 
momentum cut-off. Anyons are collective excitations which naturally have a finite size; 
however, this length does not appear in the effective Chern-Simons Lagrangian and must 
be supplemented otherwise. Possibly, it could be self-consistently determined in the exact 
solution of this theory, which is, however, not known at present. Within the mean-field 
approximation, anyons have the size given by the magnetic length O [l/^/B^j, which is the 
minimal localization of particles in the Landau levels. Therefore, we have, 



a = — 



5 



/ 5n 



5 



n 



2vrpo ' 



5 = 0(1) 



(2.37) 



where 5 is a proportionality constant. This defines the cut-off A for both the boson-based 
(9/ n = 1/ub) and fermion-based (9/n = 1 — 1/rip) anyons, 

2np A 2 27rp 



A 2 



8%n B 



A 2 



5 2 F n F 



(2.38) 



The boson-based and fermion-based ground-state energies match at the semion point 9/tt = 
1/2 for the natural choice Sp = 5b = 5; the parameter 5 = 0(1) is left free. 

We can integrate AE^ ( ^.32| ) with the respective cut-offs, use r] B = Anp /n B (respec- 
tively, t]f = 27rpo (1 + (J l n F))i an d obtain, 



E 



(2) 



E^ + AE { 



(2) 



2m 



1 + - 5 (1 + -T 

n F 2 \ n F 



n F 5 2 1 + £ 



< 



< 



-J- < I 

n B - 2 ' 



where a = tif mod 2 , and 

2 F(y) = y 2 + 2y-(y + 1) J y 2 + 2 y + log 



y + l + Jy 2 + 2y 



(2.39) 
— < 1 , 



(2.40) 



It is interesting to discuss the qualitative behavior of the ground-state energy as a func- 
tion of 9/ Ti. The quadratic correction is negative definite, as it should, and vanishes at the 
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end points 6/ir = 0,1, where the leading expressions E^ already gives the exact result. 
Near free bosons, 6/tt ~ 0, the ground-state energy is quadratic in 6/tt: this is a natural 
second-order result for bosons interacting with strength 0(9/n). Near the fermionic end, 
9/ ii = 1 — 1/uf —> 1, the oscillations 0(1 /n 2 F ) between even and odd rip values become 
subleading, and the ground-state energy approaches linearly the Fermi energy from below, 
E^ 2) ~ Eg (1 - |1 - 6/n\/25 2 ). The shape of the ground-state energy as a function of statis- 
tics has been much investigated in the quantum mechanics of a finite number N of anyons 
|2jfl . A direct comparison of these results with (|2.39|) is, however, difficult, because the 
quantum- mechanical excited states form a continuum in the large N limit. We thus find 
that the field theoretic approach gives the best result for this quantity. 



3 The U(2) non-Abelian anyon fluid 

We now consider non-relativistic matter carrying an isospin 1/2 representation of index 
r = 1,2, which interacts with an U(l) x SU(2) Chern-Simons gauge field A^ = (A^A^), 
a = 1,2, 3, with couplings = k and Ksu(2) = The Lagrangian reads 

C = z*tpo*)-^(A*) t (A*) + ^ {KdpA* - ^A^A^+^A^A, , (3.1) 
where the summation over the isospin index is implicit and the covariant derivative is 



D fl = d„ + iA fl = (d ll + iAJl + iA°— , (3.2) 



~2 

and a a are the Pauli matrices. By proceeding as in the Abelian case, we find the Hamiltonian 

H= I >x (A*) + (A*) • (3-3) 
J 2m 

The equations of motion for the Abelian part of the gauge field are again given by ( |2.3| ) and 
those of the SU(2) part are, 

_ p& = B a = €ij (diAl + l-e^A^A^ = -I p a (Gauss law), 
V 2 J J k 

F& = a A$ + d i A%-e abe AlAZ = ~e ij J? , (3.4) 
where the SU (2) isospin density and current are given by, 

H 2 

The covariant conservation law is 

{V^f = d^ a - e abc A\j^ c = . (3.6) 
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3.1 U(2) mean field approximation 



Let us look for a self-consistent approximation of the ground state which displays uniform 
densities of both matter and isospin, 

(p) = p , (Ji) = 
(P a ) = Po, (J-) = (3.7) 



Given that local gauge invariance cannot break spontaneously [|19[], a nonvanishing mean 



isospin density is not rigorously true in the exact theory. However, it is possible in 
the mean-field approximation, where local gauge invariance is reduced to the global one. 
Therefore, we shall argue that the mean field theory correctly describes the breaking of 
the U{2) global gauge symmetry down to a U{1) subgroup. Correspondingly, the isospin 
quantum number will no longer be conserved. 

We can rotate the isospin axes so that the mean isospin density become (p a ) = 5% po. 
Constant densities imply constant magnetic and iso-magnetic fields, eq. (|3]4]) and 



(Q\B a \Q) = 5 a 3 B =~6$po, (3.8) 

K 

respectively, and the corresponding fields. The zeroth-order mean-field Hamiltonian 
thus reads 

n - = | rf 2 x _L {D f)^\ D f)^ ) ( 3.9) 



where 



A (0) = ft - ,A (0) , A (0) = e/- (^1 + I a 3 ) . (3.10) 



This Hamiltonian describes two copies of Landau levels, one for each isospin orientation, 
with different values of the magnetic field B + = (B + B /2) and _B_ = (B — Bq/2), 
respectively. The number of Landau orbitals in an area AistheniVf = \B±\A/2n. Following 
the same steps of the Abelian case, we test the consistency of the mean field hypothesis 
by constructing uniform ground states for these Landau level problems. For boson-based 
anyons, we fill uniformly the pair of first Landau levels with n + isospin-up and n_ isospin- 
down particles per orbital, respectively. The equations relating densities and magnetic fields 
are 

- kB = po = 

-KBo = po = ^ |n + '~" ' 2 ' -n_ '" u o 2 1 1 . (3.11) 




12 



These equations relate the unknown quantities (n + ,n_,p ) to the external parameters 
(p Q ,K,li). In the Abelian case, the uniform filling was only possible for certain values of 
(po,n)- Here, instead, there is a one-parameter freedom, which we fix by minimizing the 
ground-state energy|. This can be expressed in terms of the known data as follows: 



E, 



(o) 



(5, 



o 



Bo\2 
2 / 



.4 



1 

2m 



4nm 
Po 

K 



•n+ + 

Po 
2k 



(B - 
Arcm 



Po , ~ 
y+Po 



-ri- 
Po 

K 



Po 
2k 



Po ~ 
y-Po 



(3.12) 



We study this expression as a function of the unknown p , with range \p \ < po/2, and 
locate its minima. Choosing for convenience k < and varying k, we find three phases for 
the non- Abelian theory: 

i) K < 0. The minimum of energy is found for p = 0; 

ii) K > and 1/|ac| > 4/|/c|. The minimum is found for po = ±po 2k/k; 

iii) k > and 1/\k\ < 4/|k|. The minimum is found for p = =tpo/2, with energy 

(o) pIa ( 1 1 



E 



po=±Po/2 



2m 



4|k| 



(3.13) 



We see that the zero-th order mean field approximation manufactures a classical potential 
with non-trivial minima, as in the standard cases of spontaneously broken symmetry (see 
fig. 1). 

The general feature of the non- Abelian problem is the presence of anyonic particles with 
both isospin charges and opposite contributions to the average iso-magnetic field Bq. If 
the two couplings have the same sign, the minimal energy ground-state configuration has 
equal populations of isospin up and down particles (n + = n_). This is the phase (i). In 
the phase (ii), the difference of populations depends on the ratio k/k. In these two cases, 
the non-Abelian mean field approximation is not self-consistent: in the first one, po = 
does not reproduce Q3.7D ; in the second case, the non- vanishing value of p is achieved in 
the singular limit n + — > oo and vanishing magnetic field in the corresponding Landau level 
(B + B 3 ) — > 0. This case is not further analyzed here. These two phases can be continuously 
connected to the theory with SU{2) gauge interaction only, by letting l/\k\ l/|?c|. In 
particular, we find that the non-Abelian mean-field approximation ( |3.7| ) is not consistent 
for Chern-Simons theories with symmetry SU(2) only, or for other semi-simple Lie algebras. 
This SU(2) invariant anyon fluid ground state is actually P and T invariant, due to the 
vanishing of the iso-magnetic field (although the fluctuations are not invariant). Similar 



models with a pair of oppositely charged anyons have been introduced |12], for explaining 



* Note that gauge invariance requires Attk to be an integer [4]. 
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1.5 



k < 




-Po/2 



-Po/4 



Po 



PO/4: 



Po/2 



Figure 1: Mean field ground-state energy (3.12) of the non-Abelian anyon fluid, plotted as 
a function of po, for one value of k in each phase. 
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the lack of P and T violations in the high-temperature superconductivity. These models 
cannot be analyzed within this mean field approximation. 

Here we shall discuss the phase (iii), where the U(l) interaction is dominating the SU(2) 
one. Actually, this phase is continuously connected with the previous Abelian model by 
letting <C In this phase, the lowest Landau level with field \B — B /2\ is empty, 

n_ = 0, and has higher energy than the other Landau level with field \B + B /2\, which is 
populated. A gap makes stable this mean-field ground state against (n + ,ri_) fluctuations. 
Note that the non-Abelian interaction lowers the ground-state energy of the pure Abelian 
theory (|2.17| ), due to the cancellation mechanism discussed above. 



3.2 U(2) current algebra 

The analysis of low-energy, quadratic fluctuations around the non-Abelian mean-field ground 
state is similar to the Abelian case in section (2.3): we must derive the non-Abelian current 
algebra and rewrite the Hamiltonian ( p.3|) in terms of currents. The non-Abelian Chern- 
Simon field can be solved completely in terms of the matter density at equal time by using 



the radial gauge [15 



x t A l (x) = . (3.14) 

These gauge conditions eliminate the commutator of non-Abelian gauge fields appearing in 
the Gauss law (p.4|), which can be solved as in the Abelian case (|2.(j|). 



A a - —— - c- .-r' T a 
~ K X ■ d 13 j ' 

4" = ^TT^^V 1 ' ( 3 - 15 ) 
k 1 + x ■ a 

where the operators (x-<9) _1 will be better defined afterwards. Actually, in (2+l)-dimensions, 
the solution of the Gauss law is possible in any axial gauge n^A^ = , at the expenses of 
breaking either rotation or translation invariance. In our problem, it is preferable to maintain 
explicit rotation invariance, because the Chern-Simons interaction is chiral. 

The non-Abelian current algebra can be obtained again by quantizing the bosonic matter 
field ^i(x) only. The commutation relations between the gauge-invariant currents (p, Jj) 
are still independent of the Chern-Simons coupling and are given by the eqs. (|2.8|) . The 
commutation relations of p a are , 

[p a (x),p(y)] = 0, 

^(x), %)] = — A (*(x - y)p a (x)) + -L e l3 yie abc p\y) (- 
J %m ox 1 imK \. 



x ■ d 
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P a (x),p"(y) 



-— z€ij x J e afec p 6 (x)5(x - y) 
imn 

i e abc p c (x) 5(x-y) . 



y ■ d 



p c (y) 



(3.16) 



We use some algebraic identities relating the matter and isospin ( 1/2 ) currents, which 
follows from the completeness of the basis of (2 x 2) isospin matrices. These are obtained 



by "gauging" the identities found in ref. |10||: 



4p a p a = p 2 



pjf = p a Ji - -e abc p\V iP ) c 



m 



(3.17) 



Therefore, we can consider the currents (p a , J») as independent variables and (p, Jf ) as 
dependent ones. In particular, we do not need the explicit form of the commutators involving 

ja 

Next, we study the normal ordering of the Hamiltonian ( |3.3| ). The analysis is similar to 
the Abelian case (eqs. (|2.20|) - (|2.21|) ) because the non-Abelian ground state has the same 
filling of the lowest Landau level. Therefore, we must add to ( p.3|) a term proportional to 
the U(2) Bogomol'nyi identity, analogous to ( |2.19| ), with coefficient a = 1. We then find: 



n 



rf 2 x 



2m 



(A*) t (A*) + ie i i(A*) t A-* 



K 



-p 



(3.18) 



This can be written in terms of the currents (p a , Jj), by using the identities (|3.17|) and 
another one for (A*) f Pi*) E 



H 



ci 2 x 



2m 



i((Ap) a (Ap) a + m 2 J 4 J 4 



[m(d l p)J j + 2 e abc £ (Ap) 6 (Ap) C j 
-p 2 -^(p a ) 2 l- (3.19) 



3.3 Quadratic fluctuations 



We expand the current algebra ( |3.16| ) and the Hamiltonian ( |3.19| ) to leading order by plugging 
in the expansions 

f>° ~ ^61 + p"(x) , J(x) ~ J(x) , (3.20) 

^4^(«T + 2T^"°) • (3 - 21) 



where the operators with hat are much smaller that their mean field values, ( see eqs. ( [2.22 ) 



( p. 24] )). The expansion of the Hamiltonian involves the covariant derivative of the isospin 
density , 

(V iP ) a ~ d iP a + ^ e^x 1 e ab3 A p\ A = 1 -— - ~ . (3.22) 
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The quadratic Hamiltonian splits into two terms (we omit the hats from now on), 

H®=Hf [p 3 ,J] + H$ A [ P \p 2 ] . (3.23) 

The first term is 



d 2 x 



2mpo 



( d iP 3 f + m 2 ( J) 2 - 'Inn, ,<),,>'. I, - p ( - + ~) (p 3 ^ 2 



(3.24) 



and is similar to the Abelian Hamiltonian ( |2.26| ). It describes local density fluctuations of 
the isospin-up particles in the lowest Landau level, without isospin flips. The second term 
in the Hamiltonian is 



-1/(2) 
n NA 



d 2 x 



2mpo 



Po 



P 



+ 



H) x 2 (Ap a )(Ap a ) 



(3.25) 



where the indices a, (3 = 1, 2 are a subset of the adjoint isospin indices. This term describes 
the non- Abelian dynamics of isospin-flip excitations. 



The approximated U(2) current algebra ( |3.16| ) of the independent variables (p a , Jj) also 
decouples into two orthogonal subalgebras of (p 3 , Jj) and (p x ,p 2 ), respectively: 



and 



P 3 (x),P 3 (y) 
P 3 (x),J(y) 



[p a (x),p Q (y)] 

p\xlp 2 (y) 



Po 



2im 







diS(x-y) 



a = 1,2 



(3.26) 



•P°x^ ^ 



(3.27) 



Therefore, the two terms of the quadratic Hamiltonian ( |3.23| ) are decoupled at the quantum 
level. The (p 3 , Jj) subalgebra is isomorphic to the U(l) current algebra (|2.23| , [2.24]) , thus 
the analysis of section 2.3 can be completely repeated for the Hamiltonian fl3.24j ). After the 
Bogoliubov transformation, one obtains a massless mode with sound velocity depending on 
both Chern-Simons coupling constants, 



V s = — 

m 



\ 



^ V \k\ 4:K , 



1 4 
< — < -; — - 



K < 



(3.28) 



Therefore, the density fluctuations which do not change isospin behaves as the Abelian ones, 
with sound velocity related to the ground-state energy (|3 
will be discussed in section five. 



). Their physical interpretation 
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In the next section, we shall find the spectrum of isospin-flip fluctuations. The (p x ,p 2 ) 
approximate current algebra can be represented by a non-relativistic bosonic field x as fol- 
lows, 



P 



X + X 



P 



M x ] - x 



x(x),x f (y)l = <*(x-y) 



(3.29) 



Their Hamiltonian can be expressed in terms of the integro-differential operator A acting 
on x, 



nj(2) 



J 2m [ 2k \2kJ 

J rf 2 x x f A 



X 



where J 



X (3.30) 

-ieijXjdi is the angular momentum and c = — po (4/ K + 1/2/c) > 0, is a posi- 
tive constant. Note that this Hamiltonian is normal-ordered, thus there is no need of the 
Bogoliubov transformation in this case; the spectrum of these excitations is given by the 
eigenvalues of the Hermitean operator A. 



4 The spectrum of non-Abelian fluctuations 



4.1 General properties and gauge invariance 



The radial gauge condition ( |3.14|) breaks translation invariance, because it selects a pre- 



ferred origin of the coordinates of the plane. This also happens in the Landau levels, due 
to the choice of the background potential Ai = ^e^x^ . Clearly, both systems have transla- 
tion invariance, which can be achieved by combining translations with compensating gauge 
transformations, the so-called magnetic translations. Non-observable non-gauge invariant 
quantities, like the eigen-functions, transform covariantly under magnetic translations, while 
physical quantities, like the energy levels, are invariant. 

The magnetic translation operators for the non-Abelian Hamiltonian ( |3.30|) are obtained 



by exploiting its close analogy with the Landau level Hamiltonian |lbfl . We introduce the 
eigen-functions ipE,i for one-particle states with energy E and angular momentum £, 

x( x ) = Hfe(x) c E ,t , (4.1) 

E,£ 

and we find that they satisfy the Schrodinger-like equation, 

A iIiej = E^ E ,i = £ D\Di + ^ ^ E ,i 

A = di- iAi = di- te^x^ ( 1 —— - -) , (4.2) 

3 2k \2 + x ■ d 2 J ' v ' 
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which is similar to the Landau-level one, were it not for the additional non-local operators 
in the covariant derivatives. In the Landau level problem, the gauge transformation which 
restores the gauge condition after translation is (B = 2), 



Xi + k, A' i -A i = d i 9, 9 = e ij x i b j , A i = ^{x? + V) 



The eigen-functions are translated and rotated 



7T; 



e- ! Vi(x + b) ~ (1 + &V + 0(6 2 ))^ L (x 
1,. i 

OXi 



while the Hamiltonian H L = —(d/dx l — iA^f /2m is covariant. This implies 

[H L ,*}] = 0. 



(4.3) 



(4.4) 



(4.5) 



In the non-Abelian problem ( (4.2|) , the compensating gauge transformation from x l Af = 
to (x i + b i )A i a = 0, with Af given by (|3.21|) , is similarly found to be 



At 
6" 



At + dr - e abc 6 b A1 

L-.xV ( s a ^ + - 

k 13 \ 3 4 + (l + x-9)(2 + a;.9) 



(4.6) 



to leading order in the fluctuating densities p a . Their approximate gauge transformations 
are found to be 



P 

X 



'3 



P 



1 _ ie .. x yEl ( I \ 

13 2k \ 2 (l + x-<9)(2 + x-<9) 



X 



(4.7) 



thus x acquires an operator-valued phase. One can show that the derivatives Dj in 
are "covariant" in the operator sense: D^e" 10 = e~ l9 Di, where 9 is given by (|4.3|) . Thus, 
in analogy with (|4.4| , [4.5|) , the operator A commutes with the following magnetic translation 
operator, 



i\i = di + ieijX 3 



i Po fl 



1 



2k \2 (l + x-d)(2 + x-d) J ' 
It is convenient to introduce the holomorphic components it, ^ , 



[A,7Ti 







7T = i(7Ti - ivr 2 ) 



d z 1 

dz + 2 ~ Z (l + x-0)(2 + a;-0) 



(4i 



(4.9) 



(hereafter, the length scale given by the non-Abelian mean field Pq/2k is put equal to 2 
and the mass m = 1). The operators 7r and tt^ satisfy the same algebra as their simpler 
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Landau-level analogues, because the non-local interaction mediated by the Chern-Simons 
field is translation invariant, 

[7T,7Tt] = l, [J,7r] = -7r, (4.10) 

These relations imply that the spectrum of A is infinitely degenerate in angular momentum, 
because 7iipE,i oc ipE,e-i or inpE/ = 0. Another analogous identity is 

2^tt = ( A - - - l) + 2 J (4.11) 



2 

which shows that tt is invertible apart from a special line in the (E, J) plane of the spectrum, 
which we shall discuss later on. Away from this line, tt^ can be used to generate the eigen- 
functions of arbitrary positive I starting from any given value, say from I = 0. 

As in the Landau-level problem, the operator A can be put into a manifest positive- 
definite form, by introducing the holomorphic components of the covariant derivatives Di 
72|), as follows, 

A = 2a f a + - - 1 , 
2 

d n l \ + d fi l \ , 

Note, however, that the eigenvalue problem is much harder, because the operators a, a) do 
not satisfy the simple harmonic oscillator algebra. The form (|4.12 ) for A allows to put a 
positive lower bound on the energy, 

which implies that the spectrum for the non-Abelian fluctuations is positive definite and has 
a mass gap. In the following discussion, we shall show that this bound is saturated and 
explain its physical origin. 



4.2 Physical interpretation of the spectrum 

In the previous section, we have shown that: 

i) the operator A looks like a non-local deformation of the Landau level problem of 
electrons in the average effective magnetic field B e g = po/2H; 

ii) the spectrum has a gap M = p (2/\n\ — I /2k), proportional to both couplings. 

Let us try to explain these results in simple terms before entering in the more technical 
analysis of the eigenvalue problem. 
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The zero-th order mean field approximation has provided us with the physical picture 
of two effective Landau level structures for up and down isospin particles, with effective 
magnetic fields pq(1/\k\ _ l/4«) and p (l/N + 1/4k), respectively. The ground state, 
corresponding to the homogeneous filling of the lowest isospin-up level, has the energy 
£o — Po(l/N — 1/4/5) /2 per particle ( p,13|) , which is given by the expectation value of 
the two-body repulsive term, 

J rf 2 x ~ (Bp + By) = S oPo A + 4£ J rf 2 x (p a f , (4.14) 

in the Hamiltonian ( |3.19| ). This repulsive interaction affects both the density fluctuations 
which are isospin diagonal (p 3 oc aia{) and isospin rotating ( p 1 + ip 2 oc ip oc aja-f ). The 
former are Bogoliubov transformed, and therefore are gapless, with sound velocity propor- 
tional to y/Eo- The latter excitations are not transformed and then acquire the gap M = 4£ 
from (14.14|) . 

In more physical terms, the isospin diagonal phonons are local fluctuations in the filling of 
the lowest Landau level, which do not feel a net magnetic field. On the other hand, the isospin 
rotating fluctuations are made of individual isospin flips, which move one electron from the 
filled up-level to an empty down-level. This can be thought of as leading to two effects: the 
hole in the filled up-level propagates as a phonon in a magnetic field, the magneto-phonon, 
which is gapful The jump of the electron to any empty down-level gives a discrete 

spectrum with steps proportional to the non-Abelian magnetic field po/2re. Therefore, we 
expect that the operator A has a discrete, Landau-like spectrum above the magneto-phonon 
gap M. 



4.3 Eigenvectors and eigenvalues 



After the separation of variables, 



^E,i(r:0)=e^y E , i (r) 



the eigenvalue problem 

2(E - A)V^(r) 



can be rewritten, for the radial part, 
1 



: ((rd r y -t) + 2£ + 2X-r 2 + 



2 + rd, 



(4.15) 



(4.16) 



where we parametrized E = A+c/2 = A+l+M. The operator in the first line of this equation 
is the Landau Hamiltonian, which would yield the spectrum A = 2k + 1, k > 0; the second 
line is the additional non-local term. It is convenient to transform (|4.16|) into a fourth-order 
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differential equation for the reduced wave-function ip = (2 + rd r )ip. By multiplying also on 
the left by (r 3 <9 r ), we obtain 

\rd r f + (r<9 r )(-r 4 + 2(A + £)r 2 - £ 2 - 4)(r<9 r ) + A£{£ + r 2 )] <p = , 

ip = (2 + rd r )cp . (4.17) 

Although the operator (2 + rd r ) has a non-trivial kernel, we shall find that it is invertible 
in the subspace of physically (normalizable) wave-functions, for which the two eigenvalue 
problems ( 4.16|) and (^4.17|) are actually equivalent. 

The analysis of the characteristic equations for the solutions of ( f4.17|) around r = and 
r = oo, leads to the following asymptotic behaviors, 







r — > oo 



■0 



3 -r 2 /2 



if, 



,-r 2 /2 
,r 2 /2 



(4.18) 



The first two behaviors for both r — > and r — ► oo in this table are found in the Landau 
problem, while the last two ones are new. Note that these four behaviors can also be obtained 
from the integro-differential form ( |4.1fl| ), by introducing two constants for the homogeneous 
solutions of the integral operators. 



Since the r — > oo asymptotics of free waves are not found in ( 4.18 ), we conclude that the 
physical solutions are square-integrable and that the spectrum is discrete. The integrable 
behaviors are ip ~ (r' £ ', r°, r 2 ) for r — » and ip ~ (e~ r2 ^ 2 , r -4 ) for r — > oo. 

Next, we analyze the action of the magnetic translation operator tt 
action of 7r generates all the integral operators 



9|). The multiple 



1 1 



r m tp(r), 



m G Z 



(4.19) 



m + rd r r m rd r 

whose action should be well-defined and unique. We first need the precise definition of these 



integral operators [{L5 



F = = /t /(Ar) • 



(4.20) 



In this equation, the constant a parametrizes the homogeneous solution of rd r F = F, which 
corresponds to the residual gauge freedom within the radial gauge fl3.14| ); we take a global 
complete gauge fixing by setting a = oo in Q4.20Q , where ip and (p vanish. This choice leads 
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to a consistent solution of the eigenvalue problem^. This gauge choice enforces the physical 
condition that matter fluctuations vanishing at infinity should not produce a gauge field 
at infinity. In this gauge, the integral operators ( [4.1 9|j are well defined on wave-functions 
with asymptotics if) ~ e _r2//2 (r — ► oo), but can be singular on if) ~ r~ 4 ; therefore, we 
neglect the latter type of solutions. We have shown that gauge invariance (the action of n) 
imposes further conditions on the physical solutions, which also ensure the invertibility of 
the relation between and if). The explicit action of the integral operators Q4.20Q on the 
basis of polynomials times the exponential asymptotics is, 



where \l/ is the confluent hypergeometric function vanishing at r = oo [24]. This is also the 
incomplete gamma function T(j3/2, r 2 /2), which is polynomial for f3 = 2,4,6,..., and has 
logarithmic terms for /3 = 0, —2, —4, .... Its expansion for r — > is, 

J_ r , e - r y 2 f + ^ + 0(r 4 )) - r(f )2^ , * 0, -2, -4, 

rd r \ i logr 2 + 0(l) , (3 = 0,-2,-4,... 

In conclusion, we accept integrable eigen-functions with asymptotics if) ~ e~ r2//2 , for 
r — > oo, and if) ~ (r°, r 2 , rl £ l) for r — > 0. The counting of free parameters is as follows: there 
are four physical independent solutions in total (three at r ~ and one at r ~ oo ), plus the 
energy, minus the four matching conditions (d/dr) n ip(r ), n — 0, 1, 2, 3, at < r < oo, and 
the normalization condition. This counting is consistent with a unique physical solution and 
a discrete spectrum. Actually, we shall find that the system of conditions is over-determined 
( —1 free parameters ), because the logarithmic solution corresponding to the degenerate 
exponents a = 0,2, at r = 0, will not be present. 



For £ — 0, the subset of physical solutions of the differential equation ( 4.17]) are also 
solutions of the following second order equation, 

[rd r f - r 4 + 2Ar 2 - 4l <p = , ip = (2 + rd r ) tp = 2 + ^ r (f) , (4.23) 

for the function 0. Indeed, the asymptotic behaviors of this equation, ~ r ±2 (r — > 0) and 
~ e ±r ' 2 / 2 (r — > 00) correspond to if) ~ r°,r 2 , and if) ~ e ±T " 2//2 , respectively. The general 
solutions of ( |4.23| ) are readily found in terms of confluent hypergeometric functions ^ and 



$: the unique regular solution at infinity is given by 



r 2 e ^/ 2 J ^(1-M;r 2 ), k = 1,2, 3, 

^ M/(1-M;r 2 ), M 1,2,3,.., 



§ Here, we do not find any analogous of the obstruction to fixing completely the axial gauge described in 
ref. [23]. 
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as a function of the energy parameter E = 2k + 2 + M. The $ solutions corresponds to the 
discrete spectrum 

E = 7 £-{k + l) + M, M = ^(^-±) , k = 1,2,3,..., . (4.25) 
2Km m \ \k\ 2k J 

(Physical units were restored in ( |4.25|) ) . Their wave-functions are polynomial, because $(1 — 
k, 3; r 2 ) truncates at the k-th term and behaves as $ ~ r 2 (r — > 0). The corresponding i[) 
functions are also polynomial, as shown by using ( (4.22| ). Therefore, these solutions are 
acceptable. 

The second type of solutions exists for the complementary continuous range of energy, 
and behaves as $ ~ r~ 2 for r — > 0. As a consequence, the corresponding ip functions are 
logarithmic for r — > 0, 

^ ~ 2r7T^Ij '° Sr2 ' ^7^ 1, 2, 3, ... (4.26) 

Although the logarithmic behavior is square-integrable, it produces a 6(r) term in the r.h.s. of 
the eigenvalue equation (|4.16| ), due to (df + ...)ip oc S(x), which cannot be accepted. Another 
reason for rejecting these solutions is that they are mapped by tt^ and n into non-integrable 
solutions ip ~ r - '^' (r — > 0) with £ = ±1, as shown in the appendix. 

In conclusion, the spectrum for £ = is discrete and given by ( |4.25| ). Let us add some 
remarks: 

i) The number of free parameters for the reduced second-order problem is equal to zero, 
because there are two physical independent solutions, plus the energy, minus two matching 
conditions and the normalization. This is correct for a discrete spectrum and show no sign 
of the £ over-determination mentioned before. 

ii) Eigenf unctions of ([4.17|) with asymptotic ip ~ r -4 (r — > oo) correspond for £ = to the 
solutions of the inhomogeneous equation [(rd r ) 2 — r 4 + 2Ar 2 — 4]0 = 1. In the appendix, 
we show that they behave as tp = 0(logr 2 )(r — > 0) and should be rejected by the same 
arguments given for eq. (|4.26| ). This is another reason for discarding this type of solutions, 
independent of (|4.19|) . 

iii) The physical solutions ( f4.25p are of the form ip — (2 + rd r )ip, with (p regular for r = 0, 
thus they have vanishing zero mode 

r roc 

d 2 xtlj = 2tt / dr«9 r (rV) = (4.27) 
J i=o J o 

as required for density fluctuations around the mean field. 
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Figure 2: Synopsis of the physical eigenstates in the angular momentum (£) versus energy 
(k) plane (The meaning of symbols (•), (o) and (x) is given in the text). 



The £ = spectrum fl4.25|) extends for £ ^ into Landau -like levels, whose eigen- 
functions can be obtained in principle by applying the 7r and ir^ magnetic translation oper- 
ators ( |4.9||4.1Cl|) to the £ = functions. For £ ^ 0, we do not have a general explicit method 
of solution and we made a case by case analysis. The properties of the £ ^ solutions 
are summarized in fig. 2 and will be briefly discussed hereafter, leaving the details to the 
appendix. The allowed values of £ are bounded from below by the line k = —£, because for 
these values 7r is not invertible, by eq. ( |4.11| ), and annihilates the physical wave functions. 

Let us first note that the condition ( |4.27 ) is trivially satisfied by £ ^ eigen-functions, 
due to their angular dependence, while it eliminates the zero mode for £ = 0. This implies 
that another line of eigenstates should exist for £ > 0, which is not connected to the previous 
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I = states. This is possible for the energy E = M, (k = — 1), because ir is not invertible at 
£ = 1, 7T0_i i = 0, where the line stops. Moreover, the bound ( 4.13 ) forbids states of lower 
energy. Indeed, such eigenstates are found by explicit analysis of the differential equation 

(EHD- 

The general properties of (i ^ 0) solutions are the following (see fig. 2). There are 
polynomial solutions, which were found by solving the four-term recursion relation with 



Mathematica These are of two types: 

i) for k > 2 even, and i > 0, 

i; = e m r^ [a + ai r 2 + ... + a k r 2k )e- r2/2 , ( points (•)), (4.28) 

ii) for k > — 1 and i > —k, t even, 

^ = e m (bo + &ir 2 + ... + b k+e/2 r 2k+i ) e"^ 2 , ( points (o)), (4.29) 



Moreover, for odd £, there are non-polynomial solutions, represented by crosses (x) in fig. 
2, which display a double power expansion: 

= e m { r W (oo + air 2 + ■••)+ fo + ^r 2 + ■ • ■} e"^ 2 , (4.30) 

which is similar to the confluent hypergeometric function These solutions can only be 
obtained by applying n or tt^ to a neighbor polynomial solution. In the appendix, we report 
a table of the polynomial eigen-functions ( |4.28|) , (|4.29| ) for the first few values of k and £, and 
give examples of the action of it. Note that no normalizable solutions are found with energy 
k = 0; thus, there is the double of the Landau-level gap between the lowest available level 
(k = —1) and all the higher ones (k = 1, 2, ...). We do not have a physical interpretation of 
this result. 

In conclusion, the complete energy spectrum is discrete and given by eq. (|4.25|) for 
k = —1, 1, 2, 3, with each level infinitely degenerate in angular momentum [I > —k). 

Let us discuss more precisely the action of the magnetic translation operators 7r and 
7T and show that it closes on these solutions, thus ensuring their translation (and gauge) 
invariance. Equations (|4.10|J4.lT| ) imply ir^iripk-k = and nir^ipk-k-i = 0: by explicit 
calculation (see the appendix), we actually find nipk-k = 0, leading to the following pattern, 

• ■ ■ ► V'fc.-A-l y i^k-k *• 1pk,-k+l ► " " " 

+ 1 - ipk,-k < — ipk,-k+i < — 

Therefore, the representation of the (tt, n') gauge algebra is not fully decomposable into 
normalizable (£ > —k) and non-normalizable (£ < —k) states: the former are mapped by 
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7r and 7T into themselves, while the latter are also mapped by 7P into normalizable ones. 
Nevertheless, the projection of the non-normalizable states to zero is consistent^]. Actually, 
the same non-decomposable representation occurs in the elementary Landau levels, because 
the action of the operators tt l and {tt l Y on the Landau wave-functions is the same as ( |4.31| ). 
Therefore, non-decomposability seems to be a rather general property of gauge invariance in 
the Hamiltonian formalism. 

We now discuss the completeness of the basis of eigen-functions we have found. For £ = 0, 
the eigen-functions ( |4.24|) are linear combinations of the polynomials {r 2n e _r2/ ' 2 , n > 0}, 
which form a complete basis for square integrable functions of (r 2 ) with vanishing zero 
mode. For £ ^ 0, the space of solutions cannot be easily defined in mathematical terms, due 
to the involved r = boundary conditions, which can be either ip ~ e^ e r° or ip ~ e tee r^. 
Therefore, the issue of completeness cannot be easily addressed for the £ ^ subspaces, 
which are, nevertheless, isomorphic to the £ = complete basis by the ir action. Note that 
the singular behavior ip ~ e *r° at r = 0, is fully acceptable, because it cancels in gauge 
invariant quantities like p a p a oc ip'ip. 

5 Spontaneous symmetry breaking and the excitations 
of the anyon fluid 

We have been describing the Abelian and non-Abelian U(2) anyon fluids, which are non- 
relativistic gauge theories of the Chern-Simons type, and we have shown the spontaneous 
breaking of the corresponding global symmetries U(l) and U(2) — > U(l), respectively. It 
is interesting to discuss the analogies and differences with the four- dimensional Yang-Mills 
theory of the Standard Model of electroweak interactions, and identify the excitations of the 
anyon fluid with the non-relativistic analogues of Goldstone and Higgs particles, if possible. 
There are two basic differences: 

i) the non-relativistic matter fields have half of the degrees of freedom of their relativistic 
counterparts, because the latter describe both particles and antiparticles; 

ii) the Chern-Simons gauge field has no propagating physical degrees of freedom and thus 
cannot lead to the Higgs phenomenon. 

Let us first recall the superfluid, which is the canonical example of spontaneous symmetry 

™ Note that our analysis does not exclude the existence of other solutions for some £ ^ isolated values. 
However, these would not be acceptable because they would not close under the action of ir. 
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breaking in non-relativistic field theory: 



Due to Bose condensation, the field acquires the ground-state expectation value = y/po, 
which breaks the U(l) global symmetry of particle number conservation. Small excitations 
around the mean field are diagonalized by the Bogoliubov transformation. This leads to a 
massless excitations with sound velocity v s oc y/gpo, controlled by the repulsive interaction; 
moreover, the ground state does not have a definite particle number, due to the Bogoliubov 
rotation. 

The usual description [111] in terms of radial (p) and phase (8) components of the field, 



\& = \/po + P e , can be easily compared to the current algebra description (|2.8| ). The 
current is represented as, 

Ji = — lm U%tf,) = , (5.2) 

m v J m 

and the approximate commutation relations ( [2.23D , (|2.24|) become, 



[Ji, p}- — W - y) — [8, p] = z5(x - y) . (5.3) 
m 

Therefore, the would-be relativistic Higgs (p) and Goldstone (8) fields are conjugate variables 
in the non-relativistic theory, the superfluid massless mode is a Goldstone particle and there 
is no non-relativistic analogue of the Higgs particle. 

The Abelian anyon fluid is very similar to the superfluid. The "microscopic" mechanism 
leading to (p) = po is not the Bose condensation - there is no macroscopic occupation of 
a single energy level, rather a macroscopic number of particles at the same energy. Never- 
theless, there is spontaneous breaking of the global U(l) symmetry, because the Bogoliubov 
transformed ground state has no definite particle number (see sect 2.3). 

One can find a closer relation with the usual superfluid by formally integrating out the 
Chern-Simons field. The resulting self-interacting matter theory can possibly reduce to ( |5.1| ) 
for small fluctuations around the saddle-point approximation. In this sense, we can consider 
the anyon superfluid as a non-relativistic example of dynamical symmetry breaking ||13|| . 
However, this is a peculiar example, where the self-interaction should have special normal- 
ordering effects, as in sec. 2.3, 2.4, which determine different effective |\1/| 4 interactions for 
boson-based (g oc 1/ac) and fermion-based (g oc 1) anyons. 

The current algebra approach is general enough to handle this non-standard mechanism of 
spontaneous symmetry breaking. By expanding again \1/ into (p, 8) components, one obtains 
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the current Jj ~ po {did — A4) /m; this is still purely longitudinal, because the Chern-Simons 
gauge field can be (locally) reduced to a pure gauge, A4 = c^a, 



Ji 



a) 



m 



(5.4) 



Therefore, the current algebra ( |2.23| ), ( |2.24j ) is the same as in the superfluid ( |5.3| ) and 
the anyon massless mode is a Goldstone particle. Note that the current, i.e. (A{ — di6), is 
the fundamental quantity of the Landau- Ginzburg theory of superconductivity f2~6|| , but it 
has there a completely different meaning, because the gauge field has transverse degrees of 
freedom. 

The U(2) non-Abelian anyon fluid can be described in similar terms. The mean field 
value (p a ) = £3 po/2 breaks the 57/(2) global symmetry to the U(l) subgroup generated by 
(1 — cr 3 )/2. Although \l/ r has vanishing ground-state value, it is still convenient to introduce 



the parametrization of the Standard Model JL3 



exp 




(5.5) 



where 6 a are the three would-be relativistic Goldstone particles and p the would-be Higgs 
one. The non-relativistic field \EV describes only two of these degrees of freedom, while the 
other two are conjugate momenta. The approximate current algebra (|3.26| ,( j3!27| ) can be 
rewritten as follows: 



Ji 



Pa 
m 



p 1 ^ po e 2 



P ^Po + P 



P 



^ 3 (x),p(y)J ~i(J(x-y) , 
^(x),0 2 (y)] ~ltf(x-y 



(5.6) 



The density p and the isospin-diagonal phase 6 3 represent the Goldstone particle as in the 
Abelian case; this excitation does not have a well-defined isospin number. The other pair of 
would-be Goldstones are conjugate variables; they do not undergo the Bogoliubov transfor- 
mation, because the isospin-down number is conserved by the remaining U(l) symmetry. The 
effective repulsive interaction (4/|/c| — 1/k) (p a ) 2 induced by the Chern-Simons field gives a 
mass gap to this excitation. 



6 Concluding remarks 

In this paper, we analyzed the U (2) Chern-Simons theory coupled to non-relativistic matter 
with isospin 1/2. We applied the mean field approximation developed in the refs. and 
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uncover a phase of the theory where the global £7(2) symmetry breaks spontaneously to the 
£7(1) one. Besides one Goldstone excitation already present in the Abelian model, we found 
a massive excitation with non-trivial non- Abelian dynamics. Therefore, the phenomenon of 
superfluidity and superconductivity, originally discovered by Laughlin , extends smoothly 
into this phase of the non-Abelian anyon fluid. 

This theory can also be consider as a toy model of the Standard Model of the elec- 
troweak interactions: in section five, we clarified the analogies and differences between the 
Chern-Simons theory and the four- dimensional Yang- Mills theory. Although there are many 
simplifying features, the non-Abelian anyon fluid is an instructive example where a non- 
Abelian gauge theory can be explicitly solved in the low-energy limit. 

Another interesting aspect of this theory is its close relation with gravity in (2 + 1) dimen- 
sion. Actually, the Chern-Simons action for the non-Abelian group ISO(2, 1) ( respectively 
50(4) ) can be rewritten as the Einstein-Hilbert action (respectively with cosmological con- 
stant) ||27|| . A suitable coupling of matter fields to gravity might lead to a theory similar to 
the non-Abelian anyon fluid: the mean-field approximation might describe a phase of semi- 
classical cosmology, where the metric is "induced" by the (dynamical) symmetry breaking 
2|. 
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A Details of the eigenvalue problem 

Examples of £ ^ eigen-functions 

The polynomial solutions ( [4.28| ), denoted by the points (•) in fig. 2, are of the form, 
^ = e «' (2 + rd r ) up , <p(r) = r l (a + a x r 2 + ■■■ + a k ^ r 2k ~ 2 ) e~ r2/2 , (A.l) 
and exist for k > 2 , even, and £ > 0. The simplest ones are listed hereafter (£ n = (£ + n)): 



k 


ao(£) Oi 


a 2 a 3 a 4 a 5 


2 


-£i 1 




4 


-£ddi 4(3^ + 8) 


-34 1 


6 


-Wddx W 3 (5£ + 12) 


-2£ 5 £ 3 (5£ + 26) 2£ 5 (5£ + 22) -5£ 5 1 



(A.2) 
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The other type of polynomial solutions, denoted by (o) in fig. 2, are of the form, 

V = e ue (2 + rd r ) <p , <p(r) = 0£ + b + 6 x r 2 + ■ ■ • + 6 fc _ 1+£/2 r 2fc+ ^ e"^ 2 
and exist for k > — 1 and £ > — fc, £ even. The simplest ones are listed hereafter: 



k 


t 


6-i 


&o 


6i 


b 2 


-1 


2 


1 








-1 


4 


4 


1 






-1 


6 


24 


8 


1 




1 


2 


-2 





1 




1 


4 


-16 


-4 





1 



(A.3) 



(A.4) 



Note that for some even, negative, values of £, these two types of solutions actually merge. 



Examples of the action of tt and tt^ 

Let us derive some of the non-polynomial solutions, denoted by (x) in fig. 2, by applying 
ii to a polynomial solution with neighbor value of t. Consider for example the lower energy 
level k = —1, and find V'-i.i" 



■0-11 OC 7T0_ 



1,2 



d z + - - z- 



1 



A2B 



l + rd r ){2 + rd r 



(2 + rd r ) ~T e 



-r 2 /2 



-r*/2 Q J26 



e ie — - e~ r ' 2/2 
rd r t 



e ie — r e 
rdr 



-rV2 = _ e ie r ( 3 r_ 
2 1 2 2 



(A.5) 



In this derivation, we used some formal properties of (rd r )~ l which follows by integration 
by parts of QL20| ), and, at the very end, its explicit form (|4.21|) . The further application of 
tt yields: 

= ( 9 * + ! - * 7TL aTToT 5t ) ~ -If r e_r2/2 = • ( A - 6 ) 
\ 2 (1 + ra r )(2 + ro r ) J r ro, 



Actually, this vanishing result can be found more easily by collecting the common denomi- 
nator 1/(1 + rd r ). Equation ( |A.6j ) verifies the closure of the gauge algebra on the physical 
solutions with energy k — —1, as indicated in the diagram fl4.31| ). Another non-trivial action 
in this diagram is given by tt^ applied to the unphysical logarithmic solution ^>_i o found in 
( |4.24| ), which has the explicit form, 



2 + rd r 1 



(log) _ ^ -T ru r S_ e -r 2 /2 

' ' rd r r 2 



rd r 



-r 2 /2 



(A.7) 
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The action of tt^ is computed as follows, 

One can similarly compute that the (£ = 0) logarithmic solutions \& in ( |4.24| ), for any 
value of k 1, 2, 3 . . ., are mapped by 7r and 7T into £ = ±1 solutions with non-integrable 
behavior ip ~ r - ^'; actually, it is sufficient to use the (r — > 0) expansion of these eigen- 
functions [2|] which is reported in ( 4.26|) . 



Solutions with r asymptotics at (r — > 00) 

The study of the asymptotic behaviors shows that these solutions are also solutions of 
the inhomogeneous reduced problem (|4.23 ), [(rd r ) 2 — r 4 + 2Ar 2 — 4]<//~ 4 ) = 1. It is sufficient 



to solve it in the case k 1, 2, 3 . . .. The inhomogeneous solutions can be obtained from the 
Green function, 

u 1 (r)u 2 (p)@(p - r) + Ul (p)u 2 (r)Q(r - p) 
a 2 (p)W(p) 

In this equation, u\ and u 2 are the two independent, homogeneous solutions vanishing at 
(r — > 0) and (r — > 00), which are given by the $ and \l/ confluent Hypergeometric functions 
in ( 4.24J) , respectively. Moreover, a 2 is the coefficient of second-order term in the differential 



equation and W is the wronskian, a 2 (p)W(p) oc p. The resulting expression for the Green 
function integrated against the source can be expanded for asymptotic values of r. For 
(r — > 00), one recover the r -4 behavior by the cancellation of the positive and negative 
exponentials of $ and For (r — > 0), one finds the behavior, 

1 



r(i-*) 



l + O r 2 +0 r 2 logr 2 e" r /2 , (A.10) 



which leads to a logarithmic behavior for ipl and a non-normalizable one for ^ ^ 
7r^^. , by the same mechanism discussed above. 



References 

[1] For a review see: F. Wilczek ed., "Fractional Statistics and Anyon Superconductivity" , 
World Scientific, Singapore, (1990); A. Lerda, "Anyons", Springer, Berlin, (1992); S. 
Forte, Rev. Mod. Phys., 64 (1992) 193. 

[2] For a review see: R. A. Prange, S. M. Girvin, "The Quantum Hall Effect", Springer 
Verlag, New York (1990). 



32 



[3] For a review see: "Strongly correlated electron matter", proc. of the Los Alamos conf. 
1993, K. S. Bedell et al. eds., Addison- Wesley, New York (1994). 

[4] For a review see: R. Jackiw and S.-Y. Pi, "Self-Dual Chern-Simons Solitons", Yukawa 
Inter. Seminar, Kyoto, July 1991, preprint CTP-2000 (1991). 

[5] A. Liguori and M. Mintchev, Comm. Math. Phys. 169 (1995) 635. 

[6] R. B. Laughlin, Phys. Rev. Lett. 60 (1988) 2677; A. L. Fetter, C. Hanna and R. B. 
Laughlin, Phys. Rev. B 39 (1989) 9676; Y.-H. Chen, F. Wilczek, E. Witten and B. 
I. Halperin, Int. J. Mod. Phys. B 3 (1989) 1001; A. L. Fetter, "Anyons and Anyon 
Superconductivity", Trieste Summer Course Low Dimensional Quantum Field Theory 
for Condensed Matter Physicists , Stanford preprint, July 1992. 

[7] R. B. Laughlin, "New Elementary Particles in Quantum Antiferromagnets and Their 
Potential Relevance to High-Temperature Superconductivity" , Summer School on Mod- 
ern Perspectives in Many-Body Theory, Stanford preprint 1994. 

[8] J. Frdhlich and U. M. Studer, Rev. Mod. Phys. 65 (1993) 733, Comm. Mat. Phys. 148 
(1992) 553. 

[9] C. A. Trugenberger, Phys. Rev. D 45 (1992) 3807; Phys. Lett. B 288 (1992) 121. 

[10] R. F. Dashen and D. H. Sharp, Phys. Rev. 165 (1968) 1857; J. Grodnik and D. H. 
Sharp, Phys. Rev. Dl (1970) 1546. 

[11] R. Feynman, "Statistical Mechanics" , Benjamin, Reading (1982); B. T.-S. Chang, J. 
Math. Phys. 12 (1971) 971. 

[12] F. Wilczek, Disassembling anyons, preprint IASSNS-HEP-91/70; D. Senechal, preprint 
LAVAL-PHY-23/92; J. L. Goity and J. Soto, Int. J. Mod. Phys. B 7 (1993) 4595; N. 
E. Mavromatos, Nucl. Phys. (Proc. Suppl.) B 33 (1993) 145. 

[13] C. Itzykson and J.-B. Zuber, "Quantum Field Theory", McGraw-Hill, New York 
(1980); 

[14] P. A. Horwaty, Phys. Rev. D 33 (1986) 407; F. Wilczek and Y.-S. Wu, Phys. Rev. 
Lett. 65 (1990) 13; G. V. Dunne, R. Jackiw, S.-Y. Pi and C. A. Trugenberger, Phys. 
Rev. D 43 (1991) 1332, (E) D 45 (1992) 3012; G. V. Dunne, "Self-dual Chern-Simons 
theories", proc. Int. Sympos. "Field theory and mathematical physics" , M. Sorak, ed., 
Korea, June 1994, preprint [hep-th/9410065| ; D. Bak, R. Jackiw and S.-Y. Pi, Phys. 
Rev. D D49 (1994) 6778; D. Bak and O. Bergman, Phys. Rev. D 51 (1995) 1994. 



33 



[15] G. Modanese and M. Toller, J. Math. Phys. 31 (1990) 452; P. Menotti and D. Seminara, 
Ann. Phys. 208 (1991) 449. 

[16] See e.g.: A. Cappelli, C. A. Trugenberger, G. R. Zemba, Nucl. Phys. B 396 (1993) 
465. 

[17] A. Cappelli, C. A. Trugenberger, G. R. Zemba, Phys. Lett. B 306 (1993) 100. 
[18] D. Z. Freedman, G. Lozano and N. Rius, Phys. Rev. D 49 (1994) 1054. 
[19] S. Elitzur, Phys. Rev. D12 (1975) 3978. 

[20] G. V. Dunne, "Edge asymptotics of planar electron densities", preprint UCONN-93-9, 
|hep-th/9311141| . 

[21] See e.g.: A. Dasnieres de Veigy and S. Ouvry, Nucl. Phys. B 388 (1992) 715; G. 
Amelino-Camelia, Phys. Rev. D 51 (1995) 2000. 

[22] M. M. Doria, R. Menikoff and D. H. Sharp, Phys. Rev. A 37 (1988) 2605. 

[23] A. Bassetto, I. Lazzizera and R. Soldati, Nucl. Phys. (1984) 319. 

[24] Erdelyi ed., "Higher transcendental functions", Krieger, Malabar (1985). 

[25] S. Wolfram, Mathematica, Addison- Wesley, New York (1991). 

[26] R. Feynman, "Statistical Mechanics" , sec. 10.12; S. Weinberg, Prog. Theor. Phys. 
Suppl. 86 (1986) 43; S. Fubini and A. Molinari, Ann. Phys. (NY) 221 (1993) 16. 

[27] E. Witten, Nucl. Phys. B 311 (1988) 46; A. Cappelli, M. Ciafaloni and P. Valtancoli, 
Nucl. Phys. B 369 (1992) 669. 

[28] For an introduction, see: A. Cappelli, Nucl. Phys. (Proc. Suppl.) 25A (1992) 54. 



34 



